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Project  Objectives 

The  project’s  aims  were  to  develop  necessary  computational  tool  to  optimize  functionally 
graded  materials.  In  particular  their  failure  characteristics  were  to  be  elucidated  under  dynamic 
conditions.  Furthermore,  new  data  processing  techniques  to  determine  key  parameters  necessary 
to  define  their  failure  conditions  were  established.  Such  techniques  to  integrate  computational 
and  experimental  means  are  indispensable  for  complex  materials  where  direct  correlations  with 
measured  data  to  material  parameters  are  often  unachievable. 


Summary  of  Results 

Recent  advances  in  materials  processing  techniques  have  enabled  the  cost-effective 
fabrication  of  inhomogeneous  materials  including  those  classified  as  functionally  graded 
materials.  Using  advanced  modeling  approaches,  the  compositions  and  material  properties  of 
graded  materials  can  be  suitably  tailored  and  designed  to  produce  performance  maximized  for 
specific  applications.  For  example,  liners  of  gun  barrels  and  lightweight  ballistic  resistant  layers 
are  potential  military  applications  of  FGM’s.  In  the  project,  a  new  criterion  for  crack  growth  of 
graded  materials  was  proposed.  Due  to  the  spatial  variations  of  material  properties,  conventional 
fracture  criteria  cannot  be  directly  applied  to  FGM’s.  This  criterion  assumes  a  certain 
relationship  between  the  separation  energy  and  the  peak  stresses  needed  to  generate  new 
surfaces.  This  was  implemented  in  a  finite  element  code  and  simulations  were  carried  out  to  test 
its  accuracy.  This  type  of  analysis  is  needed  for  the  optimization  of  FGM  performance.  It  enables 
the  characterization  of  failure  resistant  behavior  under  dynamic  loading  and  establishes  a  linkage 
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between  micro-scale  cracking/damage  and  large-scale  fracture  of  FGM  with  a  smooth  transition 
across  the  multi-scale  stages  of  failure. 

During  this  period  of  the  research  project,  several  accomplishments  were  made.  First,  the 
simulation  work  on  dynamic  FGM  fracture,  including  formulation  for  property  determination, 
was  refined  and  completed  (Nakamura  and  Wang,  2001,  Wang  and  Nakamura,  2002). 
Manuscripts  have  been  also  prepared.  Second,  experimental  investigation  was  carried  out  on  real 
FGM  with  micro-indentation,  and  our  new  procedure  based  on  inverse  analysis  was  implemented 
(Nakamura  and  Sampath,  2000,  2001;  Gu  et  al.,  2002).  The  results  confirmed  the  effectiveness 
of  this  technique.  Third,  fracture  tests  on  thick  FGM  layer  were  conducted.  Although  it  was  an 
initial  study,  a  toughening  behavior  was  observed  as  crack  propagated.  These  studies  are 
summarized  below.  More  details  are  discussed  in  the  referenced  papers. 

One  of  major  progresses  made  during  this  period  was  accurate  formulation  for  dynamic 
crack  propagation  of  FGM  using  physically  consistent  models.  Simulations  were  carried  out  with 
double  cantilever  beam  configuration  (Fig.  1).  The  proposed  criterion  assumes  the  power  law 
relationship  between  the  separation  energy  T0  and  the  peak  stress  amax  required  to  generate  new 
crack  surfaces. 


G 

max 

Here  T0*  and  Gmax*  are  the  parameters  defined  at  a  reference  point  (e.g.,  elastic  phase).  The 
exponent  a  can  be  any  real  number  that  defines  the  variation  of  Gmax  as  a  function  of  T0. 
Without  the  above  relation,  independent  measurements  of  spatial  variant  T0  and  Gmax  would  be 
extremely  complex. 
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Fig.  1.  Double  cantilever  beam  model  for  dynamic 
FGM  crack  simulations. 


Computed  results  show  highly  dependent  nature  of  crack  propagation  on  the  variation  of  Gmax 
defined  through  the  exponent  a  (Fig.  2). 

Using  the  computational  procedure  established  for  the  crack  propagation  in  graded 
materials,  a  more  complex  dynamic  failure  behavior  is  analyzed.  Here  a  FGM  model  under  high 
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velocity  rod  impact  is  considered.  Unlike  the  previous  model  with  well-defined  crack  path,  the 
impact  generates  multiple  crack  nucleations  at  arbitrary  locations  and  growths  along  arbitrary 
directions.  The  impact  load  model  consisting  of  a  protective  FGM  layer  and  a  metal  substrate 
(Fig.  3).  This  is  a  simplified  model  and  not  meant  to  simulate  actual  projectile  striking  or 
penetration  phenomena  of  real  armors.  However,  the  model  is  sufficient  to  verify  the 
explicabilities  of  the  current  failure  procedure  for  graded  materials  in  complex  boundary  value 
problems. 
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Fig.  3.  Schematic  of  protective  FGM  layer  and  substrate  subjected  to  impacting  rod. 


The  impact  generates  high  stresses  that  induce  fracturing  underneath  the  impact.  The  stress  state 
at  8ps  after  the  impact  is  shown  in  Fig.  4. 


Fig.  4.  Shades  of  stresses  at  t  =  8ps  after  impact  shown  for  pressure  and  effective  stress  Discontinuous  stresses  are 

due  to  cracking. 


In  order  to  quantify  the  effects  of  material  gradation  on  fracture  behaviors.  The  geometry, 
loading  conditions  and  material  properties  of  ceramic  and  metal  are  kept  the  same  while  different 
compositional  gradations  are  prescribed.  To  describe  the  gradation  conveniently,  the  variation  of 
material  composition  is  set  to  follow  the  power-law  relation  as, 

Vm  =  (z/d)n.  (2) 
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Here,  Vm  is  the  volume  fraction  of  metal  phase,  z  is  the  distance  measured  from  the  top  surface,  d 
is  the  thickness  of  the  FGM  layer  ( d  =  20mm),  and  n  is  the  material  gradation  parameter.  The 
metal  volume  fraction  is  always  Vm  =  0  at  z  =  0,  and  Vm  =  1  at  z  =  d.  The  crack  profiles  for  n  = 
1/3  and  3  are  shown  in  Fig.  5.  The  energy  transfer  from  the  external  work  (e.g.,  kinetic  energy  of 
rod)  to  various  internal  energy  components  including  the  dissipative  energy  (fracture  energy  plus 
plastic  flow)  are  shown  in  Fig.  4. 


Fig.  5.  Crack  profiles  at  t  =  8(is  for  n  =  1/3  and  3  gradations.  Greater  vertical  cracks  are  observed  with  larger  n. 
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Fig.  6.  Evolutions  of  various  energy  components  in  the  FGM  layer  shown  for  various  gradations  for 
n  =  1/3  and  3.  Total  dissipation  is  the  sum  of  separation  and  plastic  energies. 


In  addition,  mode  dependency  on  the  fracture  criterion  was  also  investigated.  With  larger 
Mode  II  toughness,  the  initiations  are  delayed  significantly  and  different  energy  evolutions  are 
observed.  In  all  cases,  the  energy  directly  consumed  by  crack  formations  represent  a  small  part 
of  overall  energy.  However,  the  presence  of  cracks  significantly  alters  the  other  energy 
components,  including  the  plastic  dissipation.  These  results  attest  the  importance  of  accurate 
crack  modeling  in  similar  boundary  value  problems. 

Before  actual  FGM  or  functional  materials  can  be  utilized  in  real  applications,  their 
properties  must  be  identified.  However,  due  to  the  complexity  of  the  material  system,  any 
procedures  to  quantify  the  properties  can  be  very  elaborated.  In  order  to  alleviate  such 
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difficulties,  a  new  measurement  procedure  based  on  the  inverse  analysis  was  also  introduced. 
This  procedure  uses  the  Kalman  filter  to  make  best  estimates  of  fracture  properties  of  FGMs.  The 
flow  chart  of  this  procedure  is  shown  in  Appendix 
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steel  substrate 


Fig.  7.  Micro-indentation  of  FGM  layer 


Indented  displacement  (pm) 

Fig.  8.  Close  match  between  the  simulations  and 
measurements  when  estimated  properties  are  prescribed. 


The  inverse  technique  to  determine  other  mechanical  properties  was  implemented  on  real 
specimen.  The  graded  materials  composed  of  PSZ  and  NiCrAlY  phases  were  fabricated  at  Stony 
Brook’s  Center  of  Thermal  Spray  Research.  They  were  indented  with  instrumented  micro¬ 
indentation  as  shown  in  Fig.  7.  The  measured  indented  load-displacement  record  was  used  as 
input  for  the  inverse  analysis  and  the  FGM’s  compositional  grading  as  well  as  its  effective 
properties  through  thickness  were  determined.  An  excellent  agreement  obtained  between  the 
experimental  measurements  and  the  simulated  results,  shown  in  Fig  8,  confirmed  the  accuracy  of 
this  technique  on  real  FGM.  The  cross-section  of  FGM  tested  in  shown  in  Fig.  9.  The  estimated 
grading  from  the  inverse  analysis  and  the  composition  obtained  from  a  separate  image  analysis 
are  shown  in  Fig.  10.  An  excellent  match  is  observed  in  the  results. 


Fig.  9.  Micrograph  showing  cross-section  of  thermally  Fig.  10.  Through-thickness  variations  of  PSZ  and  NiCrAlY 
sprayed  graded  coating.  volume  fractions  determined  from  the  inverse  analysis. 

NiCrAlY  variation  obtained  by  image  analysis  of  cross- 
sectional  micrograph  is  shown  with  dots 
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A  further  work  is  underway  to  determine  other  material  properties  using  a  similar 
procedure.  This  includes  fracture  parameters,  diffusion  parameters  of  composites,  elastic- 
plastic  anisotropic  materials.  The  present  results  based  on  real  specimen  offer  promising 
prospects  for  applications  of  the  inverse  analysis  to  the  above  systems.  At  Stony  Brook,  a 
state  of  art  nano-indentation  system  was  purchased  through  an  ARO  equipment  grant 
(DAAD19-01-1-0709)  with  an  equal  match  from  the  State  University  of  New  York.  The 
technique  developed  under  this  grant  will  be  utilized  for  the  new  indenter  to  probe  the 
material  characteristics  at  nano-  and  micro-levels.  In  order  to  maximize  extraction  of 
accurate  material  information  the  schematic  based  on  an  inverse  analysis  technique,  as 
illustrated  in  Fig.  11,  will  be  carried  out. 
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Fig.  11.  Schematic  of  an  inverse  evaluation  procedure  with  indentation  measurements. 


Fracture  testing  of  read  FGM  for  was  also  conducted.  In  collaboration  with 
Mechanics  of  Advanced  Materials  Laboratory  at  Stony  Brook,  thick  FGM  specimens  ( t  = 
3~5mm)  were  investigated  for  their  fracture  behaviors  as  shown  in  Fig.  12.  Since  direct 
determination  of  any  parameters  which  define  responses  of  complex/heterogeneous  materials 
is  very  difficult,  novel  data  interpretation  schemes  based  on  various  inverse  analysis 
techniques  are  being  established.  Initially  the  fracture  related  parameters  of  FGM’s  are 
estimated  from  measurable  record  and  the  degradation  of  composite  panels  is  estimated  from 
surface  strain  measurements.  As  more  complex  materials  are  considered  in  the  Army  as  well 
as  advanced  engineering  applications,  special  data  processing  schemes  to  evaluate  the 
performance  of  materials  become  more  essential. 

In  a  separate  project,  under  partially  funding  of  ARO,  stress  and  failure  of  electronic 
packaging  work  was  conducted  (Gu  and  Nakamura,  2002).  As  more  I/O  connections  are 
required  in  new  generation  of  semiconductors,  flip-chip  geometry  is  becoming  more  popular. 
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Fig.  12.  Fractured  FGM  under  three-point-bend  load.  The  load-displacement  curve  shows  the 
increasing  resistance  behavior  as  crack  propagates  into  tougher  metallic  phase 


In  this  project,  detailed  three-dimensional  finite  element  analysis  was  carried  out  for  area 
arrayed  flip-chip  packages.  The  analysis  was  used  to  quantify  the  three-dimensional  effects  on 
local  fracture  behaviors  under  thermal  load  and  to  estimate  thermal  fatigue  life  of  solder  bumps. 
Since  dimensions  of  various  components  span  more  than  three  orders  of  magnitude,  the  multi¬ 
scale  finite  element  model  approach,  were  utilized  to  analyze  the  overall  package  deformation 
and  the  local  stress  and  fracture  characteristics  near  solder  bumps.  Here  two  models  with 
different  scale,  global  and  local/cell,  were  adopted.  The  global  model  allows  identification 
critical  sites  within  the  package  where  large  deformation  occurs,  while  the  cell  model  represents 
a  small  region  of  the  package  that  contains  a  single  a  solder  bump  as  well  as  various  layers  of 
under  bump  metallurgies.  The  two-step  modeling  approach  enabled  accurate  3D  investigation  of 
interfacial  delamination  near  solder  bumps.  To  investigate  its  implication  on  fatigue  life,  the 
steady  state  cyclic  shear  was  obtained  from  repeated  thermal  loading.  The  result  shows  about 
30%  greater  shear  in  the  3D  model  and  leads  to  shorter  fatigue  life  of  package  than  the  results 
based  on  the  2D  plane  strain  analysis. 


underfill 

solder  bump 


Fig.  13.  3D  finite  element  mesh  of  the  quarter  model  and  detailed  stress  state  surrounding  the  solder  bump. 
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The  finite  element  mesh  and  internal  effective  stress  within  the  solder  bump  and  other 
nearby  components  are  shown  in  Fig.  13. 


Fig.  14.  Variation  of  /-integral  along  3D  interfacial  crack  front  between  solder  bump  and  copper  pad  is 
shown  for  the  0  =  0°  orientation  (measured  from  1  to  2).  The  solder  bump  location  at  the  comer  is 
illustrated.  The  peak  J  is  about  40%  higher  than  the  value  estimated  from  2D  plane  strain  analysis. 


In  the  analysis,  J  along  interfacial  crack  between  the  solder  and  copper  pad  was  computed.  It 
shows  strong  3D  effect  as  the  embedded  crack  tends  to  propagate  which  may  lead  to  loss  of 
electrical  connection.  This  type  of  facture  analysis  is  applicable  to  the  investigations  of 
failures  at  micro- structural  level  where  many  of  heterogeneous  materials  fail.  The  new 
procedure  developed  in  this  work  can  be  extended  many  classes  of  materials  of  the  Army’s 
interest. 
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Appendix 

The  flow  chart  of  Kalman  filter  implemented  for  the  determination  of  unknown  FGM  parameters 
is  shown  below. 


Read  in  (measured)  Afmeas  at  increment  t 


Update  xt  =  (mj,  m2) T 
xt  =  xt.1+Kt[Ars-Mxt.j )] 
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